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ABSTRACT 

The cosmic microwave background (CMB) temperature maps published by the Wilkinson 
Microwave Anisotropy Probe (WMAP) team are found to be inconsistent with the differential 
time-ordered data (TOD), from which the maps are reconstructed. The inconsistency indi- 
cates that there is a serious problem in the map making routine of the WMAP team, and it 
is necessary to reprocess the WMAP data. We develop a self-consistent software package of 
map-making and power spectrum estimation independently of the WMAP team. Our software 
passes a variety of tests. New CMB maps are then reconstructed, which are significantly dif- 
ferent from the official WMAP maps. In the new maps, the inconsistency disappeared, along 
with the hitherto unexplained high level alignment between the CMB quadrupole and octopole 
components detected in released WMAP maps. An improved CMB cross-power spectrum is 
then derived from the new maps which better agrees with that of BOOMRANG. Two impor- 
tant results are hence obtained: the CMB quadrupole drops to nearly zero, and the power in 
multiple moment range between 200 and 675 decreases on average by about 13%, causing 
the best-fit cosmological parameters to change considerably, e.g., the total matter density in- 
creases from 0.26 up to 0.32 and the dark energy density decreases from 0.74 down to 0.68. 
These new parameters match with improved accuracy those of other independent experiments. 
Our results indicate that there is still room for significant revision in the cosmological model 
parameters. 
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1 INTRODUCTION 

The CMB data from the WMAP mission is the most important basis 
of cosmology study, and the accuracies of the CMB map recovered 
from the WMAP data and its angular power spectrum are essential 
for precision cosmology. 

Recently, we find there notably exist observational effects on 
released WMAP maps. The WMAP mission measures temperature 
differences between sky points using differen tial radiom eters con- 
sisting of plus-horn and minus-horn teennett et al. 2003 ah . When 
an antenna horn points to a sky pixel, the other one will scan a ring 
in the sky with angular radius 141° to the center pixel. These mea- 
sured TOD are transformed into the full- sky temperatu re anisotropy 
map by a map-making process dHinshaw et al. 2003 al) . In released 
five-year WMAP (WMAP5) CMB maps we find significant dis- 
tortion from hot Galactic sources: the pixels in the scan ring of a 
hot pixel are systematically cooled, and strongest anti-correlations 
between temperatures of a hot pixel and it s scan-ring appear at a 
separation angle 6 ~ 141 ° iLiu &Li 2009). The above results are 
confirmed by Aurich, Lustig & Steiner (2009). Furthermore, we 
also detect the no-negligible effect of imbalance observations in 
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published WMAP5 maps jLi et al 2009): systematical dependence 
of temperature vs. observation number difference between the two 
horns of a radiometer produced by the input transmission imbal- 
ance, and significant correlation between pixel temperature and ob- 
servation number in WMAP data. These recent findings of system- 
atical error in published WMAP temperature maps push us forward 
to further check the WMAP map-making processing. 



In this work, we find a remarkable inconsistency between the 
WMAP TOD and published temperature map, which is described 
in §2. The revealed inconsistency demonstrates that there certainly 
exists a serious problem in WMAP map-making process, and it is 
worth to check the reliability of released WMAP results by repro- 
ducing CMB temperature maps from the original raw data indepen- 
dently from the WMAP team. We built a self-consistent software 
package of map-making from TOD and power spectrum estimation 
from reconstructed CMB maps. Our software successfully passes a 
variety of tests, i.e. the residual TOD test, the residual dipole com- 
ponent test, the map-making convergence test and the end-to-end 
test, which are presented in §3. With our software, improved CMB 
maps from WMAP TOD are produced and shown in §4, and the 
new angular power spectrum is derived from our new CMB maps 
and shown in §5. From the new angular power spectrum, we de- 
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termine the new best-fit cosmological parameters and present the 
results in §6. Finally, we give a brief discussion of our results in §7. 



2 INCONSISTENCY BETWEEN TOD AND RELEASED 
MAP 

Let ti denotes the temperature anisotropy at a sky pixel i. In a cer- 
tain band, the observed difference of the kth observation dk = 
t k + —t k ~, where /c + and k~ are the sky pixels pointed by the plus- 
horn and minus-horn during the observation k respectively. From 
total L observations, the differential TOD 

d = At (1) 

with A ={a(kji)} being the scan matrix. The most of elements 
a(k,i) — except for a(k,i — k + ) — 1 and a(k,i — k~) — 
— 1. The WMAP team produces the released temperature map 
i fro m the calibrated TOD d by using their map-making soft- 
ware (iHinshaw et al. 2003d : ETarosik et al. 2007h . 

Here we use a simple method to test the consistency between 
the WMAP TOD d and the released map i reconstructed from d 
by map-making. We calculate for each observation k the residual 
d k between the measured calibrated difference and the calibrated 
difference predicted by the reconstructed map, d' k — dk — (i k + — 
t k - ), to get the residual TOD 

d! = d- (£+ - i- ) . (2) 

If the temperature map is properly reconstructed, only the instru- 
ment noise should be left in d' . To check it, we produce the corre- 
lation map of the residual TOD 

to = M~ 1 A T d / , (3) 
where M = A T A is diagonally dominant jHinshaw et al. 2003ab 

M-\i,j) ~ j^8u (4) 

with Ni being the total number of observations for pixel i. If the 
temperature map published by the WMAP team is reconstructed 
correctly, the correlation map to should remain only the map- 
making error with low amplitude and no significant structured sig- 
nal on it. Because Eq. [3]is linear, we have to = M _1 A T d — 
M _1 A T (t4- — i-). For a correctly reconstructed map t , both 
M _1 A T d and M _1 A T (£ + - t_) will be equal to i , then t 
will be exactly zero, despite the inevitable numerical computation 
error, or the map making error. 

Now we check the consistency of released temperature map 
with the used WMAP TOD. In map-making, not all measured 
TOD are used by the WMAP team. Duri ng the preprocess of map- 
making, bit-coded quality flags are set dLimon et al. 20081) . A non- 
zero flag indicates that either the observation is problematic in a 
specific respect, or the beam boresight is away from one of the out 
planet s no more than 7 °, which is the antennae main beam radius 
limit (Hill et al. 2009). In both cases, the corresponding data are 
at least less optimal, and some of these data (not all) are not used 
by the WMAP team in the map-making. We follow the WMAP 
team's flagging convention^]] to get the used WMAP5 year-1 TOD. 




Figure 1. Correlation maps of residual TOD d' between the input TOD and 
the TOD predicted by the official WMAP5 year-1 temperature map of Ql- 
band, in Galactic coordinates and in units of liK. The gray area represents 
pixels without valid value. Upper panel: The data used in map making are 
those used by the WMAP team to produce the released map. Lower panel: 
The data used in map making are the measured calibrated TOD of safe- 
mode. In this mode, less data are used, see text for details of this mode. 

The upper panel of Fig.Q]shows the map to obtained by Eq.[3]from 
the used WMAP5 year-1 TOD and the released WMAP5 year-1 
temperatures of Ql-band, where visible structures along the eclip- 
tic plane and around the ecliptic poles left and the rms amplitude 
~ 265 LiK, which is much higher than the expected rms error of 
~ 0.2 /jK estimated by the WMAP t eam with flight-like s imula- 
tions for their map-making algorithm (IHinshaw et al. 2003 ah . 

More than 95% of doubtful TOD data are non-zero flagged 
for the Sun, Earth, Moon and planet avoidance, the beam boresight 
angle for one of the planets is less than 7°. However, the criterion 
really used for the maps relea sed by the WMAP team is a cut of 
only 1°.5 dLimon et al. 20081) . the revealed structure in the resid- 
ual TOD may come from ecliptic contamination by used non-zero 
flagged TOD. To test this guess, we use only the high quality data 
whose quality flags are all zero (we call this the "safe-mode") to 
produce the residual TOD and its correlation map to- The high 
quality TOD (all flags being zero) preserve more than 74% of 
all measured WMAP TOD, and the sky coverage of these data is 
more than 99.7% (neglecting the WMAP5 processing mask re- 
gion). Therefore, these data in safe-mode not only can better avoid 
the contamination from the ecliptic plane, but also are enough to 
reconstruct a reliable CMB map as well. This is also ensured by 
an end-to-end test upon the total map-making and power- spectrum 
estimation processes shown in §3.4, where they are proved to be un- 
biased and highly accurate. The to map from the safe-mode TOD, 
shown in the lower panel of Fig. [TJ still remains a structure sim- 
ilar with what from the WMAP team used TOD and a large rms 
amplitude ^291 liK. 

The above results of residual TOD test demonstrate there cer- 
tainly exists a remarkable problem in the WMAP map-making pro- 



1 (1) GENERAL FLAG - Test on bits 0, 1, 3, 4, 5. This discards data when 
the observatory is not in observing mode, the Sun is visible over the shield, 
but included data when either the Earth or Moon are visible over the shield 
(but still in the far sidelobes of the radiometer beams). (2) DA-SPECIFIC 
FLAGS - Test on bit 0, re-test on bits 1-10. Bit allows exclusion of data 



with known thermal disturbances or radiometer upsets. Bits 1-10 show that 
a planet is close to a radiometer beam. If bit 1-10 is set, compute the dis- 
tance between the indicated planet and the radiometer beam center based 
on the instantaneous pointings for all points within the frame, and discard 
only those points for which the planet lies within 1°.5 of the beam. 
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Figure 2. Correlation map of residual TOD between the input TOD and 
the TOD predicted by our new temperature map of Ql-band, in Galactic 
coordinates and in units of /iK. The gray area represents pixels without 
valid value. Upper panel: The data used in map making are those used by the 
WMAP team to produce the released map. Lower panel: The data used in 
map making are the measured calibrated TOD of safe-mode. In this mode, 
less data are used, see text for details of this mode. 

cess and using carefully tested software to improve the previous 
released WMAP results is really needed. 



3 SOFTWARE TESTING 

We write out our programs for map-making and power spectrum 
estima tion. The input of o ur map-making program is the calibrated 
TOD (iLimon et al. 2 008) downloaded from the WMAP team's 
website (ft tp://lambda.gsfc.nasa. gov/). According to the WMAP 
document (IHinshaw et al. 2003 al) . only three extra corrections are 
needed: the 1/f noise removal, the sidelobe correction, and the 
transmission imbalance correction. The sidelobe correction is ap- 
plied only to the K-band, and other two corrections are applied to 
all bands. To check the reliability of released WMAP results by 
reproducing CMB maps from the original raw data independently, 
firstly we must test our map making and power spectrum estima- 
tion softwares to make sure that they produce correct outputs. This 
includes four tests: the residual TOD test, the residual dipole com- 
ponent test, the map-making convergence test, and, finally, the end- 
to-end test checking all our data pipeline in a self-consistent man- 
ner. 



3.1 Residual TOD 

With our map-making program we produce new temperature map, 
and with the new map and Eq.[3]we calculate the residual TOD as 
already illustrated. All these are done for both the WMAP5 year-1 
Ql-band TOD and the corresponding safe-mode TOD separately, 
and their correlation maps are shown in the upper and lower panels 
of Fig. Irrespectively. In both maps, no significant structure can be 
seen and the rms amplitudes are all less than 0.15 pK, almost 2000 
times lower than what from the WMAP team's map-making prod- 
ucts and close to the rms error of ~ 0.11/xK of our map-making 
(see §3.3), indicating the inconsistent problem presented in §2 is 
prevented in our map-making algorithm. 




Figure 3. Map-making error, tout — Un (/iK), in ecliptic coordinates by 
our map-making software after 50 rounds of iterations. The input map is 
Ql-band WMAP5 temperature map without foreground removal, and the 
map-making is done with our software in safe-mode. All observations with 
either beam inside the WMAP5 processing mask are rejected in addition 
to the safe-mode data rejection. The rms of t ou t — Un in this map is ~ 
0.11/iK. The gray area represents pixels without valid value. 

It has to be noted that we repeat the map making by writing 
our ow n code but following every step given by the WMAP doc- 
ument (IHinshaw et al. 2003a) . but do not encounter their problems 
in our output maps. This demonstrates that the method of WMAP 
map-making could in principle succeed in recovering temperature 
maps, but in practice the WMAP team do not implement their own 
method correctly. 

3.2 Residual Dipole Component 

The filtered calibrated TOD can not be directly used for map- 
making, because they contain large dipole component caused by 
the joint motion of the sun and the satellite. This dipole compo- 
nent is subtracted according to the WMAP team's literature (Hin- 
shaw et al. 2008, §4). After subtracting the unwanted dipole com- 
ponent, the residual dipole component amplitude outside the KQ85 
mask in our Ql-band map is estimated to be 7.56 iiK and point to 
the Galactic coordinates (13°. 16, -11°. 78). Meanwhile, the resid- 
ual dipole component amplitude outside the KQ85 mask in the 
Ql-band WMAP5 official map is also estimated, which is 7.92 
fiK and points to the Galactic coordinates (3°. 93, -14°. 70). Both 
these dipole components are estimated with t he "remove_dipole " 
program from the standard HEALPix package jGorski et al. 2 005), 
and the difference between them is about 1.1 pK. Since the orig- 
inal dipole component amplitude is several hundreds times larger 
than the residual dipole, the neglectable difference between the two 
residual dipole components derived from the new and old CMB 
maps indicates that both our map-making and dipole component 
removal have high accuracy and their results are reliable. 

3.3 Convergence of Map-making Algorithm 

The map-making algorithm is an iterative one, thus it is impor- 
tant to test its co nvergence. An official test had been adopted by 
the WMAP team (IHinshaw et al. 2003 al) . They generated artificial 
TOD from simulated CMB map (Un) according to the scan scheme 
of the real observations, then use the artificial TOD to do map- 
making and compare the result t ou t with Un to see if their software 
converges well. The output tout is expected to be very close to the 
known input U n , and the rms of tout — tin is expected to decrease as 
the number of iterations increases, right as presented in the WMAP 
team's literature (Hinshaw et al. 2003a, Fig. 2, 50 rounds of itera- 
tions). 

We repeat this test with our map-making software for four 
individual tests: (1) Un being real WMAP temperature maps and 
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Figure 4. End-to-end test. Upper panel: the input ACDM power spectrum 
(solid line) and the final cross power spectrum (dotted line). Lower panel: 
power differences between the final cross-power spectrum and the cross 
power spectrum generated by synfast, C fina i - C synfast . 



the scan scheme being safe-mode; (2) similar to the above test but 
the scan scheme being all observations; (3) and (4) are similar to 
(1) and (2) respectively, but tin being simulated CMB maps plus 
WMAP-like foreground. Each test was run for 50 rounds of itera- 
tions. Theoretically speaking, the results of these four tests should 
be very close to each other, because the inputs and map-making 
environments do not differ much among them. The test results are 
right as expected: The maps of t ou t — tin given by all the four tests 
are similar to the WMAP team's official test (Hinshaw et al. 2003a, 
Fig. 2), and all the four rms of t ou t — tin are less than 0.2 /xK, 
which do not exceed the residual rms given by the WMAP team. 
The result of test (1) for Ql-band is presented in Fig. [3] The results 
of tests (2)-(4) are similar to Fig. [3] With such tests, we are confi- 
dent that our software is able to make reliable CMB maps from the 
WMAP TOD. 



3.4 End-to-End Test 

This is the most important test, in which the overall map-making 
and power spectrum estimation are checked in an end-to-end 
way: starting from the best-fit WMAP5 ACDM power spectrum 
(Cacdm), we use the "synfast" routine in HEALPix software 
package (available at http://healpix.jpl.nasa.gov) to generate sim- 
ulated sky maps for V and W-bands, and then use these maps 
to produce artificial TOD in safe-mode. These artificial TOD are 
used to do map-making and then the recovered maps by the map- 
making process are used to compute the final cross-power spectrum 
(Cfinai)- In the end, C 'final is compared to Cacdm to finish the 
end-to-end test. 

Such a test can definitely illustrate two things: whether the 
~ 74% observations used in map-making are enough to recover 
an unbiased CMB map, and whether our data pipeline is unbiased 
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Figure 5. Downgraded (to N side = 8) Ql-band difference maps in Galac- 
tic coordinates and in units of /iK. The gray area represents pixels without 
valid value. Upper panel: the map produced with our map-making software 
from the TOD used to make the published WMAP5 year-1 map minus the 
corresponding safe-mode map. Lower panel: The WMAP5 year-1 official 
map minus the corresponding safe-mode map. 



and accurate. If the final power spectrum is consistent with the in- 
put ACDM power spectrum, then we can say that our entire data 
processing has been tested for all aspects and been well justified. 

We give the results of this test in Fig. HI More than 99% of 
Cfinai — Cacdm is attributed to the random fluctuation intro- 
duced automatically by synfast to simulate the cosmic variance. 
To isolate and show the tiny error attributed to the overall map- 
making and power spectrum estimation processes, we derive cross 
power spectrum from the V and W-bands maps generated by syn- 
fast (C S ynfast) with the same power spectrum estimation program 
and compare Cfi na i with C sy nfast in the lower panel of Fig. |4] 
This is the true variation introduced by the overall map-making and 
power spectrum estimation processes. 

For / = 2-700, the average is only 0.0052 iik 2 and rms only 
0.51 iik 2 for Cfinai — Csynfast, the average is only —6.49 iik 2 and 
rms is 234.96 iik 2 for Cfinai — Cacdm - Therefore, we see no sys- 
tematical bias in either Cf ina i - C syn fast or Cf ina i - Cacdm, 
and the overall map-making and power spectrum estimation has 
high accuracy. In other words, our map-making and power spec- 
trum estimation has passed the end-to-end test and has been well 
justified. 

The correlation maps of residual TOD shown in §3.1 demon- 
strate that our map-making product is consistent with the input 
TOD, but the published WMAP maps contain remarkable con- 
structed error component produced during map-making. The re- 
sults of testing residual dipole component amplitude (§3.2) and 
convergence of iterative map-making algorithm (§3.3) show that 
our software has good performance. The end-to-end simulations 
(§3.4) give a prove that all our data pipeline is comprehensively 
self-consistent. We therefore expect the temperature map recovered 
from WMAP TOD by our software can provide more reliable result 
for CMB study. 

It is worth to note that, with our map-making program, we also 
reconstruct temperature maps from the TOD used to produce the re- 
leased WMAP5 map. We find that, despite a few minor differences, 
the obtained map is much more consistent with our safe-mode map 
than with the WMAP team's official map, as presented in Fig. [5] 
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Figure 6. Upper panel - Ql-band temperature maps in units of mK and in Galactic coordinates. The monopole and dipole components have been removed for 
all maps. From left to right: the official WMAP5 map; the new WMAP5 map reconstructed from the TOD used in WMAP5 map-making but by our software; 
the safe-mode map reconstructed from high quality WMAP5 TOD by our software. Lower panel - Left: the quadrupole component of the official WMAP5 
CMB map; Middle: the difference between the official WMAP5 map and the new WMAP5 map (smoothed to N side = 8); Right: the difference between the 
official WMAP5 CMB map and the safe-mode map (smoothed to N s id e = 8). 



which illustrates the safe-mode being not the reason of the major 
difference between the WMAP release and our result. 



4 TEMPERATURE MAPS 

With our map-making software we reconstruct sky temperature 
maps from WMAP TOD. In this work we discuss and solve for the 
temperature only, which is the average of all four WMAP channels: 
d = (di3 + di4 + c/23 + c?24)/4. The final temperature maps are 
derived from TOD with 80 rounds of iteratio ns. The ma p-making 
is done in N S ide — 512 resolution dGorski et al. 2005h and the 
HEALPix nested pixelization scheme only. Observations that have 
either beam inside the WMAP5 processing mask are rejected in our 
map-making, hence all pixels inside the WMAP5 processing mask 
have no values in the final temperature map. However, this does not 
affect the power spectrum estimation, because the WMAP5 pro- 
cessing mask is completely covered by the mask used for power 
spectrum estimation (the KQ85 mask or the enlarged KQ85 mask, 
which is discussed in £15. IK the pixels with no values will not be 
introduced into the power spectrum estimation. 

The upper-left plot of Fig. [6] is the WMAP5 official release 
for the Ql-band, and the upper-middle plot shows the new CMB 
temperature map of the Ql-band derived from the WMAP5 used 
TOD by our map-making software (both have the best-fit dipole 
and monopole components removed, and the plotted temperature 
ranges are -0.5 - 0.5 mK). The differences between them (official 
WMAP5 map minus our map) are computed. We downgrade the 
difference map to a low resolution map with the resolution param- 
eter N S ide = 8 and show in the lower-middle plot of Fig.[6j 

The rms temperature fluctuation on the smoothed difference 
map in the lower-middle plot of Fig. [6] is 6.6 /jK, which is much 
higher than all known map making errors. We notice that, the down- 
graded difference map in the lower-middle plot of Fig.^is almost 
the same with the WMAP5 quadruple component (lower-left of 
Fig. [6), which indicates the quadruple component in the WMAP5 
CMB map is almost completely artificial. 

We also reconstruct a safe-mode map from the WMAP5 Ql- 



band TOD of safe-mode by our map-making software and compute 
its difference to the official WMAP5 map, the results, shown in the 
upper-right and lower-right plots of Fig. [6] respectively, are very 
similar to what produced from the WMAP5 used TOD, demon- 
strate again the observed remarkable difference is mainly produced 
by the WMAP map-making process. 



5 CMB ANGULAR POWER SPECTRUM 

5.1 Estimation Method 

To produce a reliable CMB power spectrum, the foreground emis- 
sion should first be subtracted from a raw temperature map to 
make the clean CMB temperature map. Our single-year clean 
maps are produced with the foreground removal technique used 
by the WMAP team (iffinshaw et al. 20071 : iGold et al. 2009h . then 
the best-fit dipole and monopole components are subtracted before 
power spectrum estimation. After these operations, 30 V and W- 
bands single-year temperature maps are used to compute 435 cross 
power spectra with the KQ85 mask (same as WMAP team's five- 
year cross power spectra estimation scheme). For the safe-mode 
maps, due to data rejection pixels in some areas of the single-year 
maps contain no usable observation. To exclude these pixels (as 
well as pixels with too few observations) in computing the CMB 
power spectra, an extra mask that mainly lies on the ecliptic plane 
is used in addition to the KQ85 mask (which is the default mask 
for WMAP team's power spectrum estimation) to constitute an en- 
larged KQ85 mask as shown in Fig. [7] 

The released WMAP CMB power spectrum is derived from 
temperature maps via several main steps including the cross-power 
spectrum estimation, cut- sky correc tion, beam function correc- 
tion, and pixel window correction (Hinsh awetal. 2 003b). Two 
softwa re packages, HEALPix and PolSpice, a re designed to d o 
these dGorski et al. 20051 : ISzapudi et al. 200 ll : IChon et al. 2004h . 
The main difference between them is that HEALPix does not sup- 
port the cut- sky correction directly, but PolSpice does. Therefore, 
we use PolSpice in this work to produce power spectra from our 
single-year maps. HEALPix are also used with the cut- sky correc- 
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Figure 7. The enlarged KQ85 mask used in computing single year cross- 
power spectra. 




Multiple moment I 

Figure 8. Binned cross-power spectra (I = 33 - 675). Solid line: the re- 
leased WMAP5 spectrum. Dotted line: the spectrum derived by us from re- 
leased WMAP5 maps to test our power spectrum estimation processes (very 
close to the solid line and is hard to distinguish). Dashed line: the spectrum 
of new CMB maps from WMAP5 used TOD and our map-making software. 
Dashed-dotted line: the spectrum of safe-mode CMB maps from WMAP5 
safe-mode TOD and our map-making software. 



tions done by our own software, and the results are almost identical 
with PolSpice results. 
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Multiple moment L 

Figure 9. Binned cross-power spectra from raw maps before foreground 
subtraction. Solid line: from the released WMAP5 maps. Dashed line: from 
our safe-mode maps. 
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Multiple moment I 

Figure 10. Low-/ cross-power spectra (I < 30). Asterisk: the released 
WMAP5 spectrum. Circle: the spectrum of new CMB maps from WMAP5 
used TOD and our map-making software. Filled circle: the spectrum of safe- 
mode CMB maps from WMAP5 safe-mode TOD and our map-making soft- 
ware. 



5.2 Power Spectrum 

For a start, to test our power spectrum estimation processes, we 
derive a binned power spectrum from the released WMAP5 CMB 
maps with the same binning scheme as the WMAP team used. The 
obtained spectrum is nearly identical with what provided by the 
WMAP team (see the solid line and dotted line in Fig.(8j where the 
two WMAP5 power spectra obtained by the WMAP team and by us 
respectively are so close to each other that it is hard to distinguish 
them). Then we compute the binned cross power spectrum from our 
new CMB maps with WMAP5 used TOD and that with WMAP5 
safe-mode TOD respectively, the results are shown in Fig. [8] for 
high order moments (Z = 33 — 675) and Fig. [TO] for low order 
moments (I < 30) respectively. 

We also compute binned cross-power spectra from raw maps 
before foreground subtraction for WMAP5 release and for our safe- 
mode maps separately, the results are shown in Fig. [9] From Fig. [9] 
we can see that the remarkable difference between the official and 
our new spectra does not come from the foreground removal pro- 
cess. 



5.3 Low Order Moments 

Fig. [To] shows large-scale binned CMB power spectra for multiple 
moment I < 30, where asterisks mark WMAP5 release, circles and 
filled circles mask our results from new WMAP5 maps and safe- 
mode maps respectively. 

The remarkable quadruple structure in the difference map 
reminds us to see whether the large scale anomalies detected 
in the released WMAP CMB power spectrum are really cosmic 
origin or just an observation effect for some artificial reasons. 
A lack of anisotropy power on the largest angular scales has 
long been observed in WMAP CMB map s feennett et al. 20 03b: 
Iffinshaw et al. 2003bl : Iffinshaw et al. 20091) . We see from Fig. \W\ 
that the power at / = 2 further drops to nearly zero in our CMB 
maps (37.3 llK 2 from our new WMAP5 maps, 24.4 llK 2 from our 
safe-mode maps, vs. 108.7 fiK 2 from WMAP5 releasefl; mean- 
while, the power spectrum uncertainty attributed to measurement 



2 The WMAP5 power spectrum and our power spectrum here are the 
pseudo-Q based cross power spectrum. 
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Table 1: Average power spectrum (/jK 2 ) 



L 


WMAP5 


BOOMRANG 1 


This work 


220 


5716(+2.8%) 


5558(reference) 


5554(-0.1%) 


500-600 


2454 


2308 


2095 


800-950 


2427(+36%) 


17 83 (reference) 


1766(-1%) 


220-950 


2584 


2367 


2174 



1 arithmetic average of BOOMRANG98 and 03 



the WMAP release, new power spectra in high order moment re- 
gion are, i n general, better consistent with the BOOMRANG CMB 
spectrum (Ruh l et al. 20031 ; I Jones et a~2 006). Especially for the 
power spectrum near the third peak (L=800-950). The relative de- 
viation from the BOOMRANG spectrum at the first peak (Z = 220) 
is rsj 3% for WMAP5 and only ~ 0.1% for our safe-mode spec- 
trum respectively; at the third peak (/ = 800 — 950), that is ~ 36% 
for WMAP5 and only ~ 1% for our spectrum, respectively. 



errors at / = 2 is ~ 5.60 fiK 2 onl}0. The quadrupole component 
(/ = 2) drops for nearly 78% from the WMAP5 value, whereas the 
power differences between the WMAP5 and our spectra at other 
multiple moments in Fig. [10] are all < 10%. It seems unlikely that 
a nearly zero power spectrum at / = 2 can be covered by the cosmic 
variance. 

The unexplained orientation of large-scale patterns of CMB 
maps in respect to the ecliptic frame has long puzzled cosmolo- 
gists. By analyzing CMB maps from the first year WMAP data, 
Tegmark et al. (2003) found both the CMB quadrupole (I = 2) 
and octopole (/ == 3) having power along a particular spatial axis , 
and more works jde Oliveira-Costa et aL 2004; Eri ksen et al. 2 004; 



ISchwarz et al. 20041 ; Ijaffe et al. 20051) found that the preferred di 
rections of these two low-/ components are highly aligned and close 
to the ecliptic plane. Recently, significant signal spatially corre- 
lated with the ecliptic plane was detected by Diego et al. (2009) 
in published WMAP5 map and by Jiang, Lieu and Zhang (2009) 
from analyzing spectral variation. In our maps, the unexplained 
quadrupole-octopole alignment disappears; in the WMAP5 release, 
the angles between the preferred directions of CMB quadrupole and 
octopole average out at 6°. 5 for all 30 single year maps; whereas 
our single year maps give 26°. 1 average angle between the two 
components. Fig. ITTI compares the quadrupole and octopole com- 
ponents averaged over 30 V and W bands our single-year maps 
and the corresponding WMAP5 maps results directly. We can see 
that both the direction and amplitude of the quadrupole component 
have changed significantly from the WMAP5 official maps to our 
new maps. 



5.4 High Order Moments 

Fig. [8] shows the binned cross power spectra for high order mo- 
ments (/ = 33 - 675), where the solid line is the released WMAP5 
spectrum, the dashed line the spectrum of new CMB maps from 
WMAP5 used TOD and our map-making software, and the dashed- 
dotted line the spectrum of safe-mode CMB maps from WMAP5 
safe-mode TOD and our map-making software. From Fig. [8] we 
see the new power spectra being systematically lower than the re- 
leased WMAP5 power spectrum over multiple moment / = 200 - 
675. For example, by averaging all the power spectra in this range, 
both HEALPix and PolSpice give our safe-mode power spectrum 
being averagely 13% lower than WMAP5. Such a power spectrum 
decrease will cause important consequences on the best fit cosmo- 
logical parameters. 

It is interesting to see from Table 1 that, in comparison with 

3 The WMAP5 power spectrum uncertainty at I = 2 consists of the mea- 
surement errors and the cosmic variance. The later is dominating, but it has 
nothing to do with map making or power spectrum estimation uncertainties. 
Therefore, we don't have to discuss it here. 



6 COSMOLOGICAL PARAMETERS 

The peak height ratio (first peak/ second peak, which is closely re- 
lated to the baryon density) in the CMB power spectrum is ~ 2.2 
for WMAP5 release and ~ 2.6 for our spectrum, indicating that at 
least the baryon density should rise. To roughly estimate to what 
extent the cosmological parameters can be affected, a brutal force 
least- square fit search for the new best-fit parameters is done in 
the multiple m oment range from 33 to 675 by using the CMB- 
Fast package (Zaldarriaga et al. 19981 ; IZaldarriaga & Seliak 2000; 
ISeliak et al. 1996h with the ACDM model for our new WMAP5 
binned power spectrum and safe-mode spectrum separately. The 
two sets of results are close to each other. The results from the 
safe-mode spectrum are listed in Table 2. From Fig. [8] we can see 
the new power spectrum differs by no more than 20% from the 
WMAP5 power spectrum. Therefore, the WMAP5 uncertainties of 
the cosmological parameters are roughly suitable for the new cos- 
mological parameters, as adopted in Table 2. 

From Table 2 we can see that, in comparison with the WMAP5 
release, some best-fit cosmological parameters from the improved 
CMB maps are significantly different, e.g. the total matter density 
Q m = Q c + increases by ~ 25% (from 0.26 up to 0.32) and 
the dark energy density Q\ decreases by ~ 10% (from 0.74 down 
to 0.68). 



7 DISCUSSION 

In this work we observe remarkably systematic artifacts in the re- 
leased WMAP CMB map, structured signals and noises related 
with the Solar System, as shown by Fig.Q] Fig.|5]and Fig.[6j which 
turn into the almost all observed quadrupole component of the re- 
leased map. From Fig.Qj]we can see that the detected ecliptic con- 
tamination mostly comes from the map-making process of the of- 
ficial WMAP temperature map, the non-zero flagged TOD used by 
the WMAP team in their map-making are the secondary source. 
It will be very helpful if the WMAP team can thoroughly recheck 
their map-making process to find out where and how the error oc- 
curs, because it is difficult to do so for a non-WMAP team member. 

In this work, we built a self-consistent software package of 
map-making from TOD and power spectrum estimation from re- 
covered CMB map. Our software successfully passes through a va- 
riety of tests. With the tested software we reconstruct new CMB 
maps from WMAP TOD, which are significantly different with the 
official WMAP maps. Our map-making process prevent the prob- 
lem of inconsistency between input TOD and reconstructed map; 
the anomaly structure and high fluctuation observed in the residual 
TOD of official WMAP maps disappear in the residual TOD of our 
new maps. Therefore, the new CMB maps are certainly improved 
by our map-making and more reliable than the WMAP release. 

The difference between the old and new temperature maps is 
almost the same with the quadrupole component of the released 
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Figure 11. Quadrupole (upper panel) and octopole (lower panel) components averaged among 30 V and W bands single-year maps. From left to right: 
WMAP5 release, from new maps reconstructed with WMAP5 used TOD and our map-making software, from our safe-mode CMB maps. 

Table 2: The best-fit cosmological parameters 



Description 


Symbol 




Value 








WMAP5-only 1 


WMAP5+BAO+SN 1 


This work 


Hubble constant (km/s/Mpc) 


Ho 


71.91^;? 


70.1 ± 1.3 


71.0 ±2.7 


Baryon density 


n b 


0.0441 ± 0.0030 


0.0462 ±0.0015 


0.052 ± 0.0030 


Dark matter density 


n c 


0.214 ±0.027 


0.233 ±0.013 


0.270 ±0.027 


Dark energy density 


n A 


0.742 ± 0.030 


0.721 ±0.015 


0.678 ± 0.030 


Flue. Ampl. at 8h~ 1 Mpc 


erg 


0.796 ± 0.036 


0.817 ±0.026 


0.921 ±0.036 


Scalar spectral index 


n s 


963+ - 014 
u.yoo_ 015 


960+ - 014 
u.you_ 013 


0.957 ±0.015 


Reionization optical depth 


T 


0.087 ±0.017 


0.084 ±0.016 


0.109 ±0.017 


1 from Hinshaw et al. (2008) 



WMAP5 CMB map (Fig. [6]). Being consistent with this, we also 
see a structure similar to the released WMAP5 quadrupole compo- 
nent after smoothing the map of the WMAP residual TOD (Fig.[T]) 
to N S ide = 8. The correlation map of the residuals between in- 
put TOD and predicted by the released WMAP temperature map 
is also very similar to the quadrupole component observed in the 
released WMAP CMB map. Such phenomenon strongly suggests 
that, in the WMAP5 CMB temperature map, a very large portion 
of the presently known quadrupole component is actually artificial. 
As a natural consequence, in our new CMB map the / = 2 fluctu- 
ation falling near the ecliptic plane no longer exists, and then the 
quadrupole component is nearly zero. 



Despite cosmic variance, the anisotropy at the lowest harmon- 
ics as inferred from our reprocessed sky maps are likely to be 
too low to accommodate the standard cosmological parameters. A 
nearly zero quadrupole component of CMB map places a tight con- 
straint on the nature of the early universe. As well known, the sub- 
tended angle of the cosmological horizon at last scattering is only 
about 1 - 2 degree. Measured CMB power at the large angular scale 
of lowest moments should reflect the circumstance of the inflation 
period, e.g. the state of matter and the production of density fluctu- 
ations during inflation. Our result indicates that we may catch sight 
of the very early rigid universe, the primordial density perturba- 
tions should be generated in a later stage of inflation. It is hopeful to 
study the generation and revolution of density fluctuation through 



more reliable large scale CMB data and with suitable mathematical 
tool for data analysis in the future. 

Besides the ecliptic foreground contamination revealed in this 
work, we have found other two kinds of systematic distortions in 
official WMAP maps. One is the Galactic foreground contamina- 
tion. Although the WMAP team has paid huge attention to sepa- 
rate the temperature anisotropy of CMB from Galactic foreground 
emission, in released WMAP CMB maps we sti ll find notable 
distortion from hot Galactic sources (iLiu & Li 20091) : pixels 141° 
away from the hottest pixels in the Galactic plane are on average 
12-14 llK cooler than average pixels. Another one is the no - 
neglisible effect of imbalance WMAP observations dLi et al 2009h : 
significant correlation between the map temperature released by the 
WMAP team and the difference in the observation numbers that 
a pixel is pointed by the satellite's two antennae, the temperature 
distortion caused by this observational effect is as high as up to 
~ 20 liK. These two kinds of systematics are also detected in our 
new maps and should be further inspected and corrected. Therefore 
we believe that there still exist space to improve CMB maps and 
accuracies of cosmological parameters again by further improving 
map-making from WMAP TOD data and removing residual sys- 
tematic errors in recovered temperature maps. 

Except the quadrupole component, the power differences be- 
tween the WMAP5 and our low-/ spectra at other multiple mo- 
ments (Fig. [TO]) are not too large. The anomalously planar cosmic 
octopole still remains in our new WMAP5 CMB map, and also in 
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the safe-mode map just with a little lower amplitude (see Fig. [To] 
and Fig. [Til . However, from our recent analysis, we find that the 
effect of imbalance WMAP observations may contribute to some 
extent to the observed octopole component in WMAP CMB maps, 
and the real spectrum at / = 3 might be lower than observed. It is 
expected that more reliable low-/ spectrum can be derived after bet- 
ter understanding and properly correcting systematics in observa- 
tion and problems in map-making. A reliable large-scale spectrum 
will certainly help to understand the early universe. For example, 
for the primordial density perturbations, inflation theory predicts 
the index n in the scalar spectrum P(k) ~ k n being nearly unit, 
n ~ 1, over a very wide range of k by virtue of scale invariance: 
most perturbations exited the horizon with the same amplitude in 
the brief interval of rapid exponential expansion. There is now a 
possibility that a larger n at low / might be more consistent with the 
data. This in turn would mean that quantum fluctuations generated 
during inflation may no longer be a viable source of super-horizon 
perturbations without significant modification. 

The best fit cosmological parameters of this work shown in 
Table 2 are just preliminary results. It is interesting that, for nearly 
all the basic cosmological parameters (Ho, Q&, Q c , Qa, as, n s ), 
both our values and the WMAP + BAO + SN cosmological pa- 
rameter values jHinshaw et al. 2009b are simultaneou sly higher or 
lower than the WMAP5 alone JHinshaw et al. 2009b . This shows 
that the new cosmological parameters given by us (although still 
rough) match with improved accuracy those of other two different 
kinds of measurements: the distance measurements from the Type 
la supernovae (SN) and the Baryon Acoustic Oscillations (BAO) in 
the distribution of galaxies. 
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